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1. 


INTRODUCTION 


The  efforts  described  herein  are  part  of  a project  to  develop  and 
Implement  a data  processing  software  system  for  the  management  and 
analysis  of  data  received  from  the  A FGL  MESA  (Miniature  Electrostatic 
Accelerometer)  experiment  flown  aboard  three  of  the  NASA  Atmosphere 
Explorer  (AE)  series  of  satellites,  AE-C,  AE-D,  and  AE-E.  The  MESA 
data  processing  system  and  the  analytical  techniques  utilized  to  extract 
drag  accelerations  from  the  total  accelerometer  output  signal  of  each  of 
these  satellites  have  been  described  previously  (1,2). 

This  report  describes  the  development  of  an  extensive  neutral 
atmospheric  density  data  base  generated  from  the  MESA  experiments  on 
these  satellites  (taken  from  the  elliptical  phases  of  their  respective  mis- 
sions). This  data  base  represents  the  most  extensive  set  of  neutral 
density  measurements  in  existence.  Knowledge  of  atmospheric  density 
and  its  variations  is  required  for  satisfactory  low  altitude  satellite  de- 
sign, operation  and  orbit  prediction. 

Also  included  in  this  report  is  a description  of  an  empirical  model 
of  density  given  as  a function  of  altitude,  latitude,  and  geomagnetic 
activity  which  was  generated  from  this  data  base.  Further  such  studies 
by  AFGL  scientists  utilizing  this  data  base  are  in  preparation. 

1. 1 Satellite  Characteristics 

The  Atmosphere  Explorer  program  involved  three  low -altitude 
satellites  designed  to  permit  a coordinated  study  of  the  thermosphere. 

A detailed  description  of  the  AE  satellites  and  their  mission  may  be 
found  in  Reference  (3). 

AE-C,  AE-D,  and  AE-E  were  launched  into  low  perigee-high 
eccentricity  orbits.  The  on-board  propulsion  system  was  used  to  main- 
tain an  elliptical  orbit  for  approximately  one  year  for  AE-C  and  AE-E. 
AE-D  data  acquisitions  were  terminated  after  three  and  one-half  months 
when  the  spacecraft  power  system  failed.  AE-C  and  AE-E  were  put  into 
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circular  orbits  near  2f>0  km  following  their  elliptical  orbit  phases.  Data 
were  acquired  in  both  spinning  (4  rpm)  and  despun  (1  rpo)  modes. 

The  AE  data  base  as  presently  configured  and  stored  at  AFGL  in- 
cludes only  elliptical  orbit  data.  Circular  orbit  data  are  stored  at  the 
AE  computer  facility  at  NASA/GSFC. 

1.2  MESA  Experiment 

The  accelerometer  experiment  on  the  AE  satellites  provides  ex- 
tremely accurate  measurements  of  orbital  accelerations.  Three  single- 
axis Instruments  orthogonally  mounted  were  flown  on  each  AE  satellite. 

Two  sensors  were  placed  with  their  sensitive  axes  in  the  spacecraft  XY 

-■ 

plane.  These  two  instruments  were  used  to  provide  redundant  capability 
for  density  measurement,  particularly  during  the  despun  mode,  as  well 
as  to  monitor  the  orbit-adjust  propulsion  system.  The  third  sensor  was 
placed  with  its  sensitive  axis  along  the  spacecraft  Z-axis.  This  sensor 
was  used  to  monitor  spacecraft  roll  and  horizontal  winds,  particularly 
in  the  despun  mode  of  operation.  A complete  description  of  the  MESA 
experiment  may  be  found  in  Reference  (3). 

The  approach  taken  to  extract  drag  accelerations  from  the  total 
accelerometer  signal  was  provided  by  numerical  filtering  techniques. 

For  despun  orbits  atmospheric  drag  information  are  sensed  as  low  fre- 

Q A 

quency  accelerations  whose  amplitudes  vary  between  10  and  5X10 “ g's 
depending  upon  satellite  altitude.  To  measure  the  drag  information  a 
low  pass  filter  was  used.  For  spinning  orbits  the  drag  information  ap- 
pears as  a modulated  signal  at  the  spin  frequency.  A band-pass  filter 
was  designed  to  process  these  data.  A complete  description  of  the  fil- 
tering techniques  may  be  found  in  Reference  (2). 


1.3  Data  Base  Development 

The  MESA  Data  Reduction  System  (DRS)  which  has  been  developed 
and  previously  described  in  Reference  (1)  is  capable  of  extracting  raw 
telemetry  data  from  the  AE  data  base,  editing  and  temperature -correcting 
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the  telemetry  data,  extracting  atmospheric  drag  values  utilizing  digital 
filtering  techniques,  and  calculating  atmospheric  density  information. 

These  density  data  (calculated  every  quarter-second)  along  with  partic- 
ular orbital  and  attitude  Information,  are  stored  in  MESA  geophysical 
unit  (GU)  files.  In  addition,  density  values  are  interpolated  from  these 
GU  files  and  stored  in  15-second  intervals  in  unified  abstract  (UA) 
history  files  for  use  by  other  experimenters. 

The  MESA  density  data  base  was  developed  from  the  liA  history 
files  in  the  following  manner.  For  each  satellite  (AE-C,  l)  and  E), 
orbital  and  density  values  were  calculated  at  5-km  intervals  from  135  km 
to  245  km,  utilizing  Lagrangian  interpolation  techniques  and  were  stored 
in  data  files  on  the  AE  computer.  The  data  from  these  files  were  then 
shipped  to  AFGL  for  use  on  the  AFGL  CDC  5600  computers. 

At  AFGL  the  data  were  then  merged  with  Jacchia  71  and  MS1S 
model  density  values,  solar  flux,  magnetic  index  Kp,  and  geomagnetic 
latitude  and  longitude  values.  These  parameters  make  up  the  MESA 
density  data  base  for  each  AE  satellite. 

In  analyzing  the  data  from  these  data  bases  it  became  evident  that  a 
time  history  data  base  would  also  be  very  useful.  We,  therefore,  re- 
ordered the  altitude  data  base  in  chronological  format  for  each  satellite. 

This  resulted  in  a time  history  of  density  values  given  in  5-km  intervals 
for  each  orbit  of  processed  data. 

Having  distinct  data  bases  for  each  satellite  (-C,  -D,  -E),  we  next 
combined  data  from  all  three  satellites  and  generated  a composite  data 
base.  This  was  done  for  both  the  altitude  and  time  history  data  bases, 

I 

and  these  two  composite  data  bases  make  up  the  present  AE  data  base. 

The  format  of  the  AE  data  base  tapes  is  given  in  Figure  1. 

1.4  S3-1  Data  Base 


S3-1  satellite  orbital  parameters  and  accelerometer  density  values 
were  provided  to  us  in  approximately  6 -second  increments  for  each 
orbit  processed.  From  these  data  orbital  and  density  values  were 
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Data  Base  Tape  Format 
Altitude  Files 

The  data  base  is  stored  by  satellite  and  altitude  in  5-km  intervals. 
Each  tape  contains  many  files,  each  file  contains  two  types  of  records  — 
header  and  data  records.  The  tapes  are  blocked,  binary,  800  BPI. 

1.  Header  Record 

The  first  record  of  each  (altitude)  file  is  a header  record  in  the 
following  format: 

Word  Description  Format 


0.1 

word  count  (IWD=40) 

I 

0.2 

group  count  (JGP=1) 

r 

1 

satellite  ID  (AE-?  or  S3-1) 

A 

2 

experimental  name  (’MESA') 

A 

3 

satellite  altitude  - km 

I 

4 

Run  date  (YYDDD) 

I 

5 

#files*100  + file# 

I 

6-40 

blank 

R 

2.  Data  Records 

The  remaining  records  of  each  file  are  data  records  in  the  follow 

ing  format: 

Word 

Description 

Format 

0.1 

word  count  (IWD==40) 

I 

0.2 

group  count  (JGP=12) 

I 

1 

orbit  number 

I 

2 

date  - YYDDD 

I 

3 

GMT  - total  seconds 

I 

4 

GMT  - hours 

I 

5 

GMT  - minutes 

I 

6 

GMT  - sec. 

I 

7 

local  time  - hours 

I 

8 

local  time  - minutes 

I 

9 

local  time  - sec 

I 

10 

LEG  (U=upleg,  D=downleg) 

A 

11 

day /night  (D=day,  N=night) 

A 

12 

Spin/despun  (S=spin,  D=despun) 

A 

13 

geographic  latitude 

F 

14 

geographic  longitude  (+E) 

I 

15 

geomagnetic  latitude 

F 

Figure  1.  Data  Base  File  Format 
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Word 

Description 

16 

geomagnetic  longitude  (+E) 

17 

invariant  latitude 

18 

density  - gm/cc 

19 

J71  model  density 

20 

J71  model  density  (Kp  = 2) 

21 

MSIS  model  density 

22 

ratio  (measured/J71) 

23 

ratio  (measured/j71  (Kp  = 2)) 

24 

ratio  (measured/MSIS) 

25 

Ap  - daily 

26 

Ap  (6.7  hr.  lag) 

27 

Kp  (6.7  hr.  lag) 

28 

F7  (day  lag) 

29 

FB7  (81  day  average) 

30 

blank 

31 

Kp  - 0 hour  lag 

32 

Kp  - 3 hour  lag 

33 

Kp  - 6 hour  lag 

34 

Kp  - 9 hour  lag 

35 

Kp  - 12  hour  lag 

36 

average  value  of  32-35  above 

37-40 

blank 

An  EOF  follows  the  last  data  record. 

Figure  1.  Data  Base  File  Format 


Format 

I 


(Cont. ) 
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calculated  at  5 km  intervals  from  150  km  to  250  km  utilizing  Lagrangian 
interpolation  techniques  and  were  stored  in  data  files  on  the  AFGL  com- 
puter. As  with  the  AE  data  base,  the  S3-1  data  were  merged  with  J71  and 
MS  IS  model  values,  solar  flux,  magnetic  index  Kp  and  geomagnetic  lati- 
tude and  longitude  values.  This  merged  data  file  comprised  the  S3-1  data 
base. 

The  S3-1  data  base  information  was  further  merged  with  the  com- 
posite AE  data  base.  This  was  done  for  both  the  altitude  and  time  history 
files,  and  those  two  new  composite  data  bases  make  up  the  present  AE/S3 
accelerometer  density  data  base. 


2.  DATA  BASE  PROPERTIES  1 

As  mentioned  previously  the  present  density  data  base  includes 
only  elliptical  orbit  data.  Table  1 provides  orbital  characteristics  and 
dates  of  elliptical  orbit  data  acquisition  for  each  satellite.  The  altitudes 
at  which  data  were  obtained  are  shown  (shaded  areas)  in  Figure  2. 


Table  1.  Satellite  Orbital  Characteristics  and  Data  Acquisition  Periods 


Satellite 

Launch 

End  Elliptical 
Data  Acquisition 

Inclination 

Initial 

Perigee 

Initial 

Apogee 

AE-C 

Dec  73 

Nov  74 

68° 

156  km 

4000  km 

S3-1 

Oct  74 

May  75 

97° 

160  km 

4000  km 

AE-D 

Oct  75 

Jan  76 

90° 

156  km 

3800  km 

AE-E 

Nov  75 

Nov  76 

20° 

157  km 

3000  km 

2. 1 Statistical  Properties 

Density  values  calculated  at  5-km  intervals  were  stored  in  the  data 
base.  Associated  with  each  point  is  the  corresponding  J71  and  MSIS 
model  value,  ephemeris  data,  solar  flux,  and  Kp  value.  A histogram  of 
the  data  distribution  as  a function  of  geographic  latitude  for  each  satellite 

1 Portions  of  this  section  were  taken  directly  from  Reference  (4). 
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Figure  2.  Distribution  of  Density  Data 


is  given  in  Figures  3a  through  3d.  These  figures  show  the  number  of 
data  points  obtained  in  each  10°  latitude  band.  The  percentages  of  the 
total  amount  of  data  (combining  the  four  sets)  in  the  latitude  bands 
0 ± 20° , |20-50°||50-70°|,  and  1 70  - 90°  | are  47%,  25%,  18%,  and 
11%,  respectively.  The  large  amount  of  near  equatorial  data  is  due  to 
the  low  inclination  of  AE-E  and  the  precession  of  perigee  for  the  other 
satellites  as  Indicated  by  Figure  2.  The  data  above  |70°  | are  provided 
by  AE-D  and  S3-1.  Although  relatively  small  in  percentage  of  total  data, 
12,822  points  were  obtained  at  these  latitudes.  Figures  4a  through  4d 
show  the  data  distribution  as  a function  of  local  time,  given  the  number 
of  points  within  each  one  hour  interval.  The  figures  show  excellent  local 
time  coverage  provided  by  AE-C  and  -E.  The  distribution  of  data  with 
respect  to  geomagnetic  activity  is  shown  in  Figures  5a  through  5d.  These 
plots  give  the  number  of  measurements  during  which  the  three-hourly  Kp 
index  (with  a six  hour  lag)  fell  within  specific  one  unit  intervals.  Approx- 
imately 20%  of  all  the  data  occur  during  conditions  when  the  average  K 


NO.  DATA  POINTS  „ NO.  OATA  POINTS 


©•90  -70  -50  -30  -10  10  30  SO  70  90 
oeooNAPHrc  latitude  toeoj 

•'igure  3a.  Histogram  of  AE-C 
Data  Distribution  as  a Function 
of  Geographic  Latitude 
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Figure  3b.  Histogram  of  AE-D 
Data  Distribution  as  a Function 
of  Geographic  Latitude 
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Figure  3c.  Histogram  of  AE-E 
Data  Distribution  as  a Function 
of  Geographic  Latitude 


Figure  3d.  Ustogram  of  S3-1 
Data  Distribution  as  a Function 
of  Geographic  Latitude 


Figure  4a.  Histogram  of  AE-C 
Data  Distribution  as  a Function 
of  Local  Time 


Figure  4b.  Histogram  of  AE-D 
Data  Distribution  as  a Function 
of  Local  Time 


Figure  4c.  Histogram  of  AE-E  Figure  4d.  Histogram  of  S3-1 

Data  Distribution  as  a Function  Data  Distribution  as  a Function 

of  Local  Time  of  Local  Time 
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Figure  5a.  Histogram  of  AE-C 
Data  Distribution  as  a Function 
of  Geomagnetic  Activity 
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Figure  5b.  Histogram  of  AE-D 
Data  Distribution  as  a Function 
of  Geomagnetic  Activity 
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Figure  5c.  Histogram  of  AE-E 
Data  Distribution  as  a Function 
of  Geomagnetic  Activity 


Figure  5d.  Histogram  of  S3-1 
Data  Distribution  as  a Function 
of  Geomagnetic  Activity 
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ls greater  or  equal  to  4^.  In  Figures  6a  through  6d,  distribution  of  data 
with  solar  flux  is  shown.  These  plots  give  the  number  of  measurements 
during  which  the  7 value  (with  a one  day  lag)  fell  within  specific  five 
solar  flux  unit  intervals.  It  can  be  seen  that  the  data  base  was  obtained 
mainly  during  conditions  of  very  low  solar  flux.  Only  about  4%  of  the 
data  occur  on  days  with  F^  ? > 100.  An  example  of  the  data  base  is 
given  by  showing  results  obtained  at  the  altitude  of  180  km  for  each  satel- 
lite. Figures  7a  through  7d  show  these  density  data,  the  latitude  at  which 
each  data  point  was  obtained,  the  ratio  of  measured  density  to  the  J71 

model,  and  the  K index. 
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2.2  Modelling 

Reference  (4)  also  compares  these  accelerometer  density  values 
to  two  commonly  used  atmospheric  models,  Jacchia  71  and  MSIS.  Com- 
parisons were  made  as  a function  of  altitude,  latitude,  geomagnetic 
activity  and  local  time.  Frequency  distributions  of  the  data  were  de- 
scribed in  terms  of  the  mean  value  and  the  second,  third,  and  fourth 
moments  about  the  mean. 

Neither  model  gives  a completely  adequate  description  of  the  ob- 
served atmospheric  variability  as  a function  of  all  these  atmospheric 
parameters.  Hence  we  developed  a multiple  linear  regression  analysis 
program  to  attempt  to  generate  an  empirical  model  of  atmospheric  den- 
sity variations  in  terms  of  latitude,  longitude,  local  time,  geomagnetic 
activity,  solar  flux,  season  and  semi-annual  effects. 

In  our  multiple  linear  regression  expression  the  dependent  vari- 
able, atmospheric  density,  is  approximated  by  a linear  combination  of  a 
group  of  independent  variables,  that  is: 

n 

0 +XNVei  * 

i=l 

The  regression  technique  finds  the  coefficients,  a.,  i = l n which  mini- 

2 

mize  £ e^  that  is,  which  best  describe  the  influence  of  x.  onp. 
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The  empirical  model  is  a function  of  latitude,  Kp,  solar  flux,  an- 
nual and  semi-annual  effects,  and  it  was  run  at  selected  satellite  altitudes 
160,  180,  200,  220, and  240  km.  The  following  table  displays  statistical 
data  related  to  our  regression  analysis  of  AK-C  data.  Regression  coeffi- 
cients (Sj)  are  given  at  each  altitude.  The  multiple  regression  coefficient 
(R),  which  is  a measure  of  the  amount  of  variation  in  the  data  which  is  de- 
scribed by  the  model,  is  also  given. 

Figure  8 displays  AE-C  measured  density  and  empirical  model 
values  displayed  as  a function  of  GMT  at  200  km. 

The  present  program  forms  a basis  for  the  extraction  and  analysis 
of  finer  density  variations  in  the  AE-C  data.  Additionally  the  program 
should  be  expanded  to  include  the  entire  data  base  (AE-C,  AE-D,  AE-E 
and  S3-1)  and  may  be  easily  modified  if  more  independent  variable  terms 
are  required. 


3.  OTHER  FUNCTIONS 

3. 1 Orbital  Bias  Determination 

The  accelerometer  is  electrostatically  suspended  and  force  re- 
balanced. Ground  calibration  requires  sufficient  suspension  force  to 
operate  in  the  earth's  lg  field. 

An  important  characteristic  of  accelerometer  performance  is  its 
null  bias,  or  output  reading  in  the  presence  of  no  acceleration  input. 
Instrument  bias  results  from  cross  coupling  of  the  suspension  axis  forces 
into  the  sensitive  axis.  To  accurately  determine  atmospheric  density 
values,  the  bias  acceleration  must  be  eliminated  from  the  total  sensor 
output.  Preflight  bias  values  are  necessarily  calculated  using  a lg  sus- 
pension force.  For  each  constrainment  range,  the  sensor  is  rotated  until 
zero  output  pulse  rate  is  obtained.  The  sensitive  axis  is  then  rotated 
exactly  180° . Bias  is  computed  from 


Bias  (ng)  = ^>U^se  r^te  {PPS)  x constrainment  range  scale  factor  , 
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Tabic  2.  Regression  Analysis  Results 


ai 

a! 

ai 

ai 

ai 

i 

Fi 

160  km 

180  km 

200  km 

220  km 

240  km 

Kp(12) 


10. 7 

6 sin  a tD 

7 cos  atD 

8 sin  2a 

9 cos  2a  tj 

10  sin  OJt^ 

11  cos  u>tT 


12  sin  2cotr 


13  | cos  2totL 

Number  of 
points 

Mean  density 
<P) 

Standard  devi- 
ation of  p 

Multiple  cor- 
relation 
coefficient 


1.13E-12 

6.46E-14 

9.70E-16 

2.17E-15 

2.71E-15 

-4.38E-15 

-4.29E-14 

-2.56E-14 

-1.15E-13 

-1.01E-13 

-1.27E-14 

-1.18E-14 

-2.38E-14 

8.84E-15 


3.81E-13 

2.56E-14 

7.05E-15 

6.39E-16 

1.55E-15 

-1.31E-15 

-1.53E-14 

-3.16E-15 

-4.94E-14 

-3.66E-14 

-5.03E-15 

-5.66E-15 

2.58E-15 

6.73E-15 


3.80E-14 

1.38E-14 

3.81E-15 

2.08E-16 

1.04E-15 

6.51E-16 

-2.02E-15 

4.43E-15 

-2.12E-14 

-1.68E-14 

-6.74E-15 

-8.04E-15 

3.60E-15 

3.06E-15 


-3.56E-14 

7.58E-15 

2.72E-15 

1.04E-16 

7.79E-16 

6.56E-16 

-8.35E-16 

5.29E-15 

-1.01E-14 

-1.00E-14 

-7.64E-15 

-6.73E-15 

3.31E-15 

1.71E-15 


1.30E-12  5.35E-13  2.50E-13  1.26E-13 
0.23E-12  0.96E-13  0.50E-13  0.32E-13 


-5.03E-14 
4.18E-15 
1.72E-15 
6.37E-17 
4.58E-16 
6.58E-16 
-3. 03E-16 
2.87E-15 
-4.22E-15 
-6. 62E-15 
-4.47E-15 
-6.09E-15 
4.60E-16 
3.34E-15 

422 

6. 80E-14 
1.77E-14 


lu 

where  p = aQ  + ^ ^ a^Fj 
i=l 

Kp(12)  = 12  averaBe  of  Kp  index  previous  to  density  measurement 

Kpg  = Kp  6 hrs  before  density  measurement 

LAT  = geographic  latitude 

Fjq  rj  - 10.7  cm  solar  flux  of  previous  day 

FjQ  7 =81  day  centered  average  of  F1Q  7 

a = 2ff/365.2422  tD  = days  since  1 Jan.  1974 

co  = 2ff/24  t^  = local  time  in  hours  of  day 


AE-C  ZOO  KM 

empirical  model  oewgitt  (om/co 
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where  bias  is  given  in  terms  of  the  value  calculated  using  a lg  suspension 
level.  Extrapolation  to  low-g  suspension  levels  is  made  by  assuming 
bias  is  directly  proportional  to  suspension  level.  Bias  temperature  de- 
pendence was  determined  from  data  calculated  at  30°  F,  73°  F and  120°  F. 
Values  of  bias  at  73°  F and  the  bias  temperature  coefficient  are  summa- 
rized in  Table  3.  These  values  are  in  units  of  ^g/cross-axis  g;  hence, 


Table  3.  Bias  and  Bias  Temperature  Coefficient  Ground 
Calibration  Data  for  Atmosphere  Explorer 
Accelerometers 


AE-C 

AE-D 

AE-E 

B 

bt 

B 

bt 

B 

bt 

XY 

YX 

ZZ 

+54.3 

+36.5 

+81 

+0.05 

-0.24 

+0.28 

+103 

+35 

+110 

+0.20 

-0.99 

-0.73 

+112 
+ 52.5 
+ 93.8 

0 

-0.26 

-0.46 

B - bias  in  units  Mg/g  , where  g = suspension  force  in 

g units.  s s 

B = bias  temperature  coefficient  in  units  \ig/g  °F. 

* s 


.3 

for  a 10  suspensi  n voltage  the  maximum  expected  orbital  bias  for  each 
instrument  is  three  orders  of  magnitude  below  the  ground  deduced  values. 
The  ground  calibration  of  bias  (B)  at  any  temperature  T is  determined 
from  the  bias  temperature  coefficient  (BT)  using 

Bias  (T)  = Bias  (73°  F)  + BT(T-73°F)  . 

Postlaunch,  the  null  bias  errors  are  considerably  reduced  by  using 
a suspension  force  several  orders  of  magnitude  lower  consistent  with  the 
expected  orbital  accelerations.  However,  bias  values  in  a low-g  orbit 
environment  may  be  different  from  those  deduced  from  ground  calibration. 
Bias  was  calculated  as  the  average  value  of  48  data  points  obtained  above 
360  km  on  both  the  upleg  and  downleg  portion  of  elliptical  orbits.  Data 


are  presented  for  the  XY  instrument  on  AE-C,  and  the  YX  instrument  on 
AE-D  and  AE-E. 
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Table  4 summarizes  the  orbital  data  from  the  three  instruments 
studied.  This  table  gives  bias  values  at  73°  and  the  standard  deviations 
(in  parenthesis),  the  bias  temperature  coefficients,  and  the  factor  by 
which  orbital  bias  was  reduced  as  a result  of  lowering  the  sensor  sus- 

_3 

pension  voltage  to  10  g. 


Table  4.  Orbital  Bias  and  Bias  Temperature  Coefficient  Results 


AE-C 

AE-D 

AE-E 

Bias 

(Standard  Deviation) 

1.270  fig 
(0.133  jug) 

-1.130  ng 
(0.268  ng) 

-4.632  ug 
(0.054  Ug) 

Bias  Temperature 
Coefficient 

-0.0020  /ig/°  F 

-0.0335  tlg/°  F 

-0.0205  /ig/°F 

Reduction  Factor 
from  lg  value 

42.7 

31.0 

11.3 

No.  data  points 

400 

300 

334 

The  long  term  stability  of  the  orbital  bias  data  and  its  temperature 
dependence,  indicate  that  the  above  determined  values  can  be  used  for 
reduction  of  density  data  when  routine  bias  measurements  are  not  possible. 
These  results  show  that  bias  must  be  determined  in  orbit  rather  than  from 
ground  calibration  for  accurate  drag  measurements  in  despun  orbits. 

For  a more  complete  description  of  AE  accelerometer  orbital  bias 
determination  see  Reference  (5). 

3.2  Correlative  Measurements 

Density  data  base  development  was  previously  described  in  Section 
1.3.  In  order  to  correlate  long  term  density  variations  with  other  experi- 
ments on  AE,  we  were  required  to  create  similar  data  bases  from  two 
AE  experiments  - Open-Source  Neutral-Mass  Spectrometer  (OSS)  and 
Neutral -Atmosphere  Composition  Experiment  (NACE).  This  was  carried 
out  utilizing  the  same  methods  as  described  in  Section  1.3  using  the  OSS 
and  NACE  unified  abstract  history  files  on  the  AE  computer. 
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Composition  data  were  calculated  at  10-km  intervals  from  140  km 
to  200  km  for  AE-C  utilizing  interpolation  techniques.  These  data  were 
stored  in  data  files  on  the  AE  computer  and  later  shipped  to  AFGL  for 
use  on  the  AFGL  computer  system. 

At  AFGL  OSS  composition  data  were  me: ged  with  MESA  density 
data,  total  mass  density  values  were  calculated  and  comparisons  to  MESA 
values  were  made.  Figures  9 to  11  illustrate  some  of  the  comparisons 
which  were  made.  Further  such  studies  are  contemplated. 


3.3  Circular  Orbit  Processing 

The  objectives  of  the  AE  missions  was  to  study  phenomena  in  the 
atmosphere  at  altitudes  above  120  km.  This  was  to  be  accomplished  in 
two  orbital  phases:  elliptic  orbit  and  circular  orbit.  After  the  initial 
elliptic  orbit  phase  was  completed,  AE-C  and  AE-E  were  circularized 
for  the  second  phase  of  their  mission. 

For  circular  orbit  spinning  data  our  processing  techniques  to  cal- 
culate density  (Reference  (2))  were  for  the  most  part  unchanged.  For 
despun  circular  orbits,  since  instrument  bias  values  are  required  and 
the  satellite  altitude  was  below  the  "non -drag  regions"  used  in  the  elliptic 
orbit  phase,  an  alternate  method  was  required  to  determine  bias  values. 
We  utilized  the  techniques  described  in  Section  3.1  above  and  Reference 
(5)  to  calculate  bias  temperature  coefficients  for  each  satellite.  Circu- 
lar orbit  temperature  data  was  then  used  to  calculate  instrument  bias 
values,  and  filtering  techniques  were  used  to  determine  density  values. 
Circular  orbit  density  data  are  being  calculated  in  a production  environ- 
ment and  stored  in  the  AE  unified  abstract  files. 

One  method  of  analyzing  circular  orbit  data  is  to  normalize  calcu- 
lated density  values  to  any  altitude  utilizing  model  atmosphere  theory. 
This  is  done  by: 


normalized  density  = 


measured  density  x model  density 
(at  the  normalized  altitude) 
model  density 


DENSITYC  G-M/CC  )X10-13  0 DENS  I TY/CCX 1 09  N2  DENS  I T Y/CCX 1 09 


DENSITY! GK/CC  ) X 1 0“ 1 3 0 DENS  I TY/CCX 1 09  N2  DENSITY/CCX 


RE-C  200  KM  0SS/MESR  COMPARISON 
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Figure  12  is  an  example  of  this  type  of  result.  Normalized  density 
data  are  given  as  a function  of  local  time  and  satellite  altitude.  In  this 
example  the  Jacchia  71  model  atmosphere  was  used  to  normalize  ellip- 
tical orbit  data  to  an  altitude  of  180  km,  but  any  other  reference  model 
atmosphere  may  be  used. 


Figure  12.  AE-E  Normalized  Density  Data  as  a Function  of  Local  Time 
and  Latitude. 

3.4  EUVS  Experiment  Support 

Utilizing  our  knowledge  (gained  in  supporting  the  MESA  accelero- 
meter experiment)  of  the  AE  data  system,  the  AE  Sigma-9  computer  sys- 
tem, the  AE  Data  Management  Facility  (DMF)  software,  and  the  general 
concepts  and  philosophy  of  the  AE  Aeronomy  Team  of  scientists  we  were 
required  to  support  the  efforts  of  the  solar  ultraviolet  experiment  (EUVS) 
scientists  at  AFGL.  We  utilized  the  remote  terminal  dedicated  to  EUVS 
to  provide  interactive  on-line  support.  This  included  maintenance  and 
execution  of  interactive  and/or  batch  programs;  copying  various  output 

' t 

files  to  the  remote  printers;  promotion  and  demotion  of  telemetry, 
orbit/attitude,  geophysical  unit  and  unified  abstract  files;  and  general 
account  and  file  maintenance  to  insure  smooth  execution  of  standard  EUVS 
processing.  In  addition,  manual  plotting  and  hand  calculations  were  re- 
quired and  they  were  done  by  RDP,  Inc.  personnel. 
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4.  CONCLUSIONS 


An  extensive  data  base  has  been  developed  utilizing  accelerometer 
measurements  obtained  with  four  low  altitude  satellites.  In  addition  to 
density  values  the  data  base  incorporates  appropriate  satellite  orbital  in- 
formation, solar  and  geophysical  parameters  and  atmospheric  model 
values.  It  is  anticipated  that  this  data  base  will  be  useful  for  correlative 
studies  of  aeronomical  problems  requiring  knowledge  of  the  neutral  at- 
mosphere, as  well  as  an  improved  model  of  the  lower  thermosphere. 
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